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Abstract 

We performed the deepest search for an X-ray emission iine between 0.5 and 7 keV from 
non-baryonic dark matter with the Suzaku XIS. Dark matter associated with the Miiky Way 
gaiaxy was seiected as the target to obtain the best signai-to-noise ratio. From the Suzaku 
archive, we seiected 187 data sets of biank sky regions which were dominated by the X-ray 
diffuse background. The data sets were from 2005 to 2013. Instrumentai responses were 
adjusted by muitipie caiibration data sets of the Crab Nebuia. We aiso improved the technique 
of subtracting iines of instrumentai origin. These energy spectra were weii described by X-ray 
emission due to charge exchange around the Soiar System, hot piasma in and around the 
Miiky Way and superposition of extra-gaiactic point sources. A signai of a narrow emission 
iine was searched for, and the significance of detection was evaiuated in consideration of the 
biind search method (the Look-eisewhere Effect). Our resuits exhibited no significant detection 
of an emission iine feature from dark matter. The 3a upper iimit for the emission iine intensity 
between 1 and 7 keV was ^ 10“^ photons cm~^ s~^ sr“^, or ^ 5 x 10“"^ photons cm~^ s“^ 
sr“i per Mq pc“^, assuming a dark matter distribution with the Gaiactic rotation curve. The 
parameters of steriie neutrinos as candidates of dark matter were aiso constrained. 
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1 Introduction 

The existence of “dark matter” (DM) in the Universe has been demonstrated by modem astrophysical 
and cosmological observations, such as rotation curves of spiral galaxies, masses of intra-cluster 
gas, gravitational leasing of clusters of galaxies, and observational data of the Cosmic Microwave 
Background anisotropy and Large Scale Structure. It is considered to constitute about a quarter of 
the total energy density in the present Universe and to occupy more than 80 % of the total mass 
density. The current “standard cosmology”, the A-Cold Dark Matter (ACDM) model, assumes that 
dark energy and DM play prominent roles in gravitational effects in the formation of structure in the 
Universe. 

Since nucleosynthesis models of the early Universe place limits on the fraction of baryons, 
DM cannot be baryonic matter, which forms a part of the Standard Model (SM) of particle physics. 
The nature of DM is thus an issue of fundamental significance. Several candidates such as “WIMPs” 
(weakly interacting massive particles), “super-WIMPs” and “sterile neutrinos” postulated, and either 
direct or indirect DM searches have been conducted to find them. However, none of them have been 
successful at present. 

Searching for photon emission from the decay or annihilation of DM particles through astro- 
physical observations is a promising approach to the discovery of DM, and the X-ray region is one 
of the possible windows. In particular, sterile neutrinos as DM candidates attract a lot of attention 
(Abazajian et al. 2001a; Dolgov & Hansen 2002; Boyarsky et al. 2012). They could be generated suf¬ 
ficiently in the early Universe through given mechanisms, for example the non-resonant or resonant 
production (Dodelson & Widrow 1994; Shi & Fuller 1999) and decay of scalar field (e.g. Kusenko 
2006; Shaposhnikov & Tkachev 2006; Petraki & Kusenko 2008). The relic sterile neutrino abundance 
from scattering-induced conversion of active neutrinos was first analytically estimated by Dodelson & 
Widrow (1994) and able to account for all of DM. The model containing sterile neutrinos is strongly 
motivated by the neutrino flavor oscillation which is supported by the atmospheric neutrino evidence 
of the Super-Kamiokande (Fukuda et al. 1998), the baryon asymmetry of the Universe and other cu¬ 
rious things beyond the SM (Asaka & Shaposhnikov 2005; Asaka et al. 2005). The flavor oscillation 
between sterile neutrinos and active neutrinos or radiative decay is predicted although its mixing an- 
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gle may be very small. Assoeiated with the deeay, a photon with the energy E = msC^/2 are emitted 
(mg is the mass of sterile neutrinos, c is the speed of light). A keV-mass sterile neutrino is a Warm 
DM (WDM) eandidate (Abazajian et al. 2001b). It resolves several ineonsisteneies between the pre- 
dietions of the CDM model and the observational results sueh as the shape and smoothness of DM 
halos (Goerdt et al. 2006; Gilmore et al. 2007; Wyse & Gilmore 2008; Lovell et al. 2014). Sinee the 
keV-mass sterile neutrinos should deeay and produee a keV X-ray photon, a seareh for this radiative 
deeay emission line in the X-ray range is meaningful. 

Radiatively deeaying DM sueh as sterile neutrinos have been vigorously searehed for by X-ray 
observatories. Boyarsky et al. (2007), Boyarsky et al. (2014a), Horiuehi et al. (2014) and Malyshev 
et al. (2014) searehed DM in the Loeal Group, i.e. the Milky Way galaxy, its satellite dwarf galaxies 
and M31, and obtained the tightest restrietion on the DM line intensity. In 2014, several papers 
reported a possible X-ray emission line around 3.5 keV. In Bulbul et al. (2014), the first report of 
this line, XMM-Newton observational speetra of 73 elusters of galaxies were used to seareh for a 
non-baryonie emission line. The signifieanee of this deteetion was as high as 4.3a. Independently, 
Boyarsky et al. (2014a) found this line in the Perseus eluster and M31 with the 4.4a signifieanee. 
Other groups (Tamura et al. 2015; Carlson et al. 2015), however, eould not eonfirm this deteetion with 
the same instruments or the same targets and gave the upper limits. At the eenter of the Milky Way 
galaxy, Boyarsky et al. (2014b) reported to deteet the 3.54 keV line by XMM-Newton, while Riemer- 
Sprensen (2014) and Jeltema & Profumo (2015) obtained upper limits by Chandra and by XMM- 
Newton, respeetively. Due to the eomplexity of the instrument responses and possible astrophysieal 
emission (e.g. Kxvill around 3.5 keV), a seareh for DM requires eareful studies. 

In this paper, we seleeted the best observational target for the DM seareh in the keV energy 
range with X-ray CCD instruments and searehed for a non-baryonie X-ray signature to provide valu¬ 
able eonstraints on the parameter spaee of DM. All error ranges in our text and tables eorrespond to 
90 % eonfidenee levels, and vertieal error bars in our figures indieate la. Throughout this paper, we 
assume the following values for eosmologieal parameters: the energy density parameters Dm = 0.27, 
Da = 0.73 and the Hubble eonstant ho = 0.7. We utilized Ftools in HEAsoft version 6.15 and XSPEC 
version 12.8.D for this analysis. 


‘ http://heasarc.gsfc.nasa.gov/lheasoft/ 
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2 Strategy of the analysis 

2.1 Selection of targets and instruments 

Selecting suitable targets and instruments is important in searching for DM in the keV energy range. 
Gravitational sources with high DM column density such as galaxies and clusters of galaxies are 
plausible candidate targets. Since we are located within the DM distribution of the Milky Way, we 
also have the potential to detect its signal over the whole sky, even in the blank sky. Such targets, 
however, contain high temperature plasma and so are sources of X-rays as background emission. 

The DM column density in the Perseus cluster or a rich cluster of galaxies was estimated as 
~ 800 Mq pc“^ (Bulbul et al. 2014; Tamura et al. 2015). That of M31 is ~ 600 (200 - 1000) Mq pc“^ 
(see Boyarsky et al. 2014a and references therein). We estimated the distribution of the DM within 
the Milky Way by the NFW profile with parameters from Sofue (2012), as po = 1-06 x 10“^ Mq pc“^ 
and h = 12.5 kpc, and calculated the column density as a function of an angle from the Galactic center 
(GC). It is ~ 200 Mq pc“^ towards the GC and decreases to 40 Mq pc“^ at the anti-Center. 

The sensitivity of the line detection is limited by the background emission. We simulated 
mock Suzaku spectra for the Perseus cluster, M31, and a typical blank sky observation with the same 
field of view (FoV) and exposure time. The blank sky is occupied by the X-ray diffuse Background 
(XDB) which consists of plasma emission in the Milky Way and a superposition of extragalactic point 
sources. Typically, the surface brightness of the Perseus or M31 are ~ 10000 or ~ 100 times brighter, 
and the line detection limits for them are ~ 100 or ~ 10 times higher than that for the blank sky. The 
meaning of “typical” is that there are much line emission from the associated plasma, which hinder 
the DM detection at various energies even if they are modeled well. Potential sensitivity for the DM 
detection is evaluated by the DM column density divided by the line detection limit. The sensitivity 
for M31 and the blank sky are almost at the same level, which is ~ 10 times better than that for the 
Perseus cluster. 

We then evaluated the sensitivity based on the actual observation cases by X-ray observa¬ 
tories, Suzaku, XMM-Newton and Chandra, with realistic instrumental responses and Non-X-ray 
Background (NXB), as well as available exposure time and sky coverage of archival data. The to¬ 
tal exposure time for archival data of the Perseus cluster is larger than that of M31, but could not 
compensate the sensitivity limited by the cluster emission as background. In the archives, there are 
hundreds of observational data sets of blank sky regions or maskable faint compact sources which 
are dominated by the XDB and used to search for DM associated with the Milky Way. The total 
photon statistics for these XDB data is larger than that for the other targets including M31. Due to 
the low surface brightness of the XDB, however, the NXB would also limit the sensitivity. The NXB 


of XMM-Newton (MOS, PN) and Chandra (ACIS-I, ACIS-S) for corresponding sky coverage are 
comparable to surface brightness levels between that of M31 and XDB observations. Thus the line 
sensitivity for the XDB observation by XMM-Newton and Chandra are limited by their NXB. On 
the other hand, the X-ray Telescope and X-ray Imaging Spectrometer of Suzaku (XRT-XIS; Mitsuda 
et al. 2007; Serlemitsos et al. 2007; Koyama et al. 2007) combined to give the lowest and most-stable 
NXB owing to combination of the low-Earth orbit and the instrumental design. Using Suzaku XIS 
observational data of the XDB is therefore the most suitable for this study. 

We note that we gave up to use the GC regions because there are number of bright point 
sources and diffuse structure which are difficult to be modeled and subtracted. In this paper, we used 
cleaned event files of the FI-CCDs (XISO, 2 and 3) and the BI-CCD (XISl) of Suzaku. 

2.2 Archival data selection 

We collected approximately-pure XDB data in the Suzaku XIS archive from 2005 to 2013 that satis¬ 
fied the requirements shown below: 

1. No bright sources (i.e. original observational purpose are blank sky fields study) or maskable faint 
compact sources are in the XIS FoV. 

2. Galactic latitudes |6| > 20° to avoid the X-ray emission peculiar to the Galactic disk (Masui et al. 
2009). 

3. Separate from the region occupied by the North Polar Spur and other local X-ray sources. 

4. The XIS was operated in the normal clocking mode and the 3 x 3 or 5 x 5 editing mode. 

Ultimately, 187 Suzaku XIS observational data sets were selected as shown in Fig. 1. We collected 
more XDB data than has been observed previously by the Suzaku XIS. Since many complex unre- 
solvable emission appears in the low energy range (<0.7 keV) and XDB photons are not adequately 
available (NXB photons dominate) in the high energy range (> 5 keV), we used the 0.7 - 7 keV and 
the 0.5 - 5 keV ranges for the FI-CCDs (XISO, 2 and 3) and the BI-CCD (XIS I), respectively. Signals 
from anomalous (hot and flickering) pixels in the XISs were screened out with the Ftool cleans is. 

3 Data reduction and reproduction of instrumental responses 

3.1 Data reduction and point source removal 

We first conducted an imaging analysis over all of the data in order to reject resolvable X-ray point 
sources contaminating the XDB spectra. The procedure was the following: 

1. XIS images of the 0.5-7 keV range were extracted. 
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Fig.l. 187 regions used to search for a keV signature of DM from Suzaku archival data. These are plotted on the all sky map with the Galactic coordinate 

system centered at the Galactic anti-center. 


2. Point sources in the XIS FoVs were deteeted and rejeeted by using a wavelet funetion of similar 
size to the point spread funetion of the XRT-XIS using wavdetect from the Chandra Interaetive 
Analysis of Observations, version 4.6. 

3. Point sourees whose flux was larger than 1.0 x erg em“^ s“^ in the 0.5 - 7 keV range were 
removed with eireular regions eentered at their positions. The radius of the eireular regions > 1.5' 
were determined sueh that these regions ineluded > 90 % of source photons. 

Seeond, in order to remove X-ray emission from the Earth’s atmosphere and from the Solar wind, and 

to reduee the NXB, we also seleeted good time intervals by applying the following eriteria, as was 

done in Sekiya et al. (2014b): 

1. The elevation angle from the bright/dark Earth limb was ehosen as > 20°/5° to avoid fluoreseent 
emission lines from the Earth’s atmosphere. 

2. Time intervals exeluding the South Atlantie Anomaly passage 

3. The Cut Off Rigidity (COR2) > 8 GV c~^ to reduee the high-energy-eharged-partiele baekground 
due to the low Earth’s magnetie field (Tawa et al. 2008). 

4. Time periods when the proton flux in the Solar wind fell below the typieal threshold, 4.0 x 10® 
em“^ s“^, to lower effeets of the Geoeoronal-Solar wind eharge exehange (Eujimoto et al. 2007). 
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The proton flux was observed with monitoring satellites: ACE/SWEPAM^ and WIND/SWE^. 

After this data sereening, the total exposure time beeomes 31.5 Msee ~ 1 year, and we obtained 
~ 2 X 10® eounts of photons in the 0.5 - 7 keV range. This is the largest data for the XDB with the 
Suzaku XIS. 

3.2 Nominal instrumental responses and NXB reproduction 

We then reproduced instrumental responses at the time of each observation. The energy redistribution 
matrix files (RMEs) including quantum efficiency and energy responses (energy scale and resolution) 
of the XISs were generated by the Etools xisrmfgen. The ancillary response files (AREs) involving 
angular responses and effective area of the XRT-XIS modules were produced by the Suzaku cali¬ 
bration database and Monte-Carlo ray-tracing simulations: the Etools xissimarfgen (Ishisaki et al. 
2007). In xissimarfgen, we assumed a uniform sky centered at each observational coordinate with 
radius of 20 arcmin as an X-ray emitting region for the simulation. 

The XDB spectra are contaminated by the NXB which originates from charged particles, elec¬ 
trical noise, and scattered and fluorescent X-ray emission from instrumental elements. The NXB con¬ 
tributions in the given spectra was also estimated from accumulated night-Earth observations around 
150 days centered at the day of each observation with the Etool xisnxbgen (Tawa et al. 2008). 

3.3 Response correction with the Crab Nebula calibration 

The nominal responses deserved above are based on the calibration data and the ray-tracing simula¬ 
tion. However, a certain level of deviation exists in the response reproduction in the 1.5 - 3.5 keV due 
to difficulty of dealing with the complicated multiple absorption edges in the XRT-XIS effective area. 
In order to check this reproducibility, we analyzed the stacked spectrum of the Crab Nebula (hereafter 
simply Crab), the data with the lowest statistical uncertainty of all Suzaku XIS data whose total count 
of photons in the 0.5 - 7 keV range is ~ 4 x 10^. The Crab emission is well modeled by a featureless 
spectrum of synchrotron emission. We assumed the emission model for the observed spectrum to be 
a power-law function absorbed by the interstellar medium of the Milky Way. We fitted the Crab spec¬ 
trum with the model multiplied by the response and found the residuals of up to ~ 10 % especially 
in the 1.5 - 3.5 keV range. The deviations were regarded as being due to a mis-reproduction of the 
response, and these were used to correct the responses. 


- http://www.srl.caltech.edu/ACE/ASC/level2/lvl2DATA.SWEPAM.html 
^ http://web.mit.edu/space/www/wind.data.html 
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3.4 Reproduction of instrumental line contributions 

The NXB including the instrumental emission lines have been estimated and subtracted with night- 
Earth observational data as described in Tawa et al. (2008), with a typical error of 3 %. There are 
known instrument lines of Al-Ka, Si-Ka, Au-Ma, Mn-Ka and Mn-K/3. We found that subtleties 
resulting from the subtraction of these lines can be mitigated by spectral fitting with a five-Gaussian 
model. 

4 Spectral analysis of the XDB 

Once the 187 sets of the XDB spectra, the RMFs, the ARFs and the NXB data were obtained, we 
then stacked the XDB spectra using the Ftool mathpha. This was done to increase the number of 
photons and so to reduce the statistical uncertainty. With consideration for both the changes in the 
instrumental-calibration and the difference among the four detectors that make up the XIS, we divided 
the entire period (2005 - 2013) into 8 periods to sort by instrumental calibration conditions (e.g. SCI 
setting) as shown in Table 1 and stacked the XDB data in each of the shorter periods. After correcting 
for systematic deviations, the total-photon-count-weighted average (in the 0.5 - 7 keV range) of 
RMFs X ARFs (responses) was made by the Ftool addrmf. The exposure-time-weighted average of 
the NXB data were also stacked. In total, 25 stacked data sets over eight periods (1 period x XISO, 
1, 2, 3 -I- 7 periods x XISO, 2, 3) were made. 

We then carried out a spectral analysis of the 25 stacked XDB spectra after subtraction of the 
NXB, using the corrected responses. We fitted the spectra in the 0.5 - 7 keV range with the typical 
XDB emission model of the following components: 

(1) Heliospheric Solar Wind Charge Exchange (H-SWCX) and Focal Hot Bubble (FHB) 

(2) Milky Way hot plasma Halo (MWH) 

(3) Cosmic X-ray Background (CXB) 

(4) Unresolved High Temperature Plasma (UHTP) 

(1): the H-SWCX is due to interaction between the Solar wind and neutral atoms in interplanetary 
space (Cox 1998; Cravens 2000; Koutroumpa et al. 2006; Yoshitake et al. 2013). The FHB is consid¬ 
ered to be the high temperature plasma with the temperature of kT 0.1 keV (T ~ 10® K) and the 
density of riH ~ 0.005 cm“^ embedded in a ~ 100 pc cavity of the cold interstellar medium in which 
the Solar System resides (Cox & Anderson 1982; McCammon & Sanders 1990). The H-SWCX and 
FHB blend is empirically represented by an unabsorbed optically-thin collisionally-ionized (CIE) 
thermal plasma emission model for the current energy range and the responses of the X-ray CCD 
cameras. (2): the MWH is considered to be hot plasma bound in the Milky Way with the temperature 
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Table 1. Stacked data properties. 


Period 

Date 

Total exposure* 

Total count^ 

Average 

2005-2006 

2005/10/01 - 2006/09/30 

3.2 

205071 

0.029 


(SCI operation started for all XlSs from October, 2006.) 


2006-2007 

2006/10/01 -2007/08/31 

4.2 

261725 

0.035 

2007-2008 

2007/09/01 -2008/08/31 

3.4 

212512 

0.029 

2008-2009 

2008/09/01 -2009/08/31 

4.4 

284447 

0.030 

2009-2010 

2009/09/01 -2010/05/31 

4.0 

242187 

0.030 

2010-2011 

2010/06/01 -2011/05/31 

4.5 

271709 

0.029 


(Injection charge increased to 6 keV for XISl on June 1, 2011.) 


2011-2012 

2011/06/01 -2012/05/31 

2.3 

250229 

0.027 

2012-2013 

2012/06/01 -2013/07/01 

3.5 

220777 

0.034 


Notes. 

* Exposure time (XISO+1+2+3) in unit of Msec after data screening, 
t Total photon count in the 0.5 - 7 keV range. 

The exposure-time-weighted average of the neutral hydrogen column density in unit of 10^^ cm“^ derived from the LAB 
Galactic H i Survey. 


of kT = 0.2 - 0.4 keV (Yao et al. 2009; Hagihara et al. 2011; Sakai et al. 2014). It is described by 
an absorbed optically-thin CIE thermal plasma emission model. (3): the CXB believed to be super¬ 
position of numerous faint extragalactic point sources such as active galactic nuclei and represented 
by an absorbed power-law model with its photon index F ~ 1.4 (Mushotzky et al. 2000; Kushino 
et al. 2002; Moretti et al. 2003). (4): the UHTP emission sometimes appears in the XDB spectra 
with strong emission of Fe-L complex and Ne-K lines (Yoshino et al. 2009; Sekiya et al. 2014a). It 
is described by an absorbed optically-thin thermal CIE plasma emission model with the temperature 
of kT = 0.4 - 1.2 keV. For the optically-thin thermal CIE emission model, we used APEC (Smith 
et al. 2001; Foster et al. 2012). The temperature of the APEC model for the H-SWCX + EHB was 
fixed to kT = 0.1 keV (Yoshino et al. 2009). The element abundances of the three APEC models were 
set to the Solar-neighbor values (Anders & Grevesse 1989) except for the Ne and Mg abundances of 
the UHTP The redshifts of the three APEC models were fixed to zero. In the power-law emission 
models for the CXB (powerlaw in XSPEC), their photon indices were permitted to vary around 1.4. 
In Suzaku XIS observational data, an OI fluorescent line from the Earth’s exosphere was sometimes 
found especially after 2011, at the time of Solar maximum, despite the fact that its contamination 
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was mostly removed by applying the elevation angle criteria as described above (Sekiya et al. 2014b). 
When we found a discrepancy in the spectral fitting result below 0.6 keV, we added a Gaussian with 
a centroid of 0.525 keV for Ol (gaussian in XSPEC). To summarize, the following model for the 
spectral fitting was adopted: [“APECi” + “Galactic absorption” x (“APEC 2 ” + “APEC 3 ” + “CXB”) 
+ “O l”] where “Galactic absorption” was for a photoelectric absorption by the interstellar medium of 
the Milky Way (phabs in XSPEC) which were able to estimate from accurate observational data of the 
neutral hydrogen column density Ah (Kalberla et al. 2005), APECi, APEC 2 and APEC 3 correspond 
to the H-SWCX + EHB, the MWH and the UHTP, respectively. 

We fitted the stacked XDB spectra with the model for each period. The results of the XDB 
spectral fitting and the best-fit parameters of the XDB model are shown in Eig. 2 and Table 2. As the 
goodness of fit, x^/dof (dof) was 1.05 (3693), and the null hypothesis probability was 2.5 %. The 
deviations between the data and model were about 1.6, 2.9 and 4.0 % in the root-mean-square at 0.5 
- 2, 2 - 5 and 5-7 keV, respectively. 

5 Search for a keV signature of DM 

We simultaneously fitted the stacked XDB spectra with the XDB model with a Gaussian emission 
line component corresponding to a DM signature. The line intensity of all the spectra were linked 
together and allowed to vary freely. The intrinsic line width of the Gaussian was assumed to be 0 eV 
because the detection of velocity dispersion due to Galactic rotation or intrinsic to WDM less than 
several 100 km s“^ are impossible with the current instrument. Its center energy was fixed and was 
allowed to vary in 25 eV steps over the 0.5 - 7 keV range. There were 261 trial lines, and these were 
not completely independent because the line profiles were broadened by the XIS energy response 
and these steps were smaller than the energy resolution in the 0.5 - 7 keV range. We determined 
the Gaussian normalization, i.e. the line intensity and its la statistical error range. The results are 
summarized in Eig. 3 (a) by blue crosses. The unit of the line intensity is shown by “Eine Unit” (EU) 
defined as photons cm~^ s“^ sr“^. The 3a upper limit on the DM line intensity was calculated as a 
sum of the best-fit line intensity and its 3cr statistical error range, or the only 3a statistical error range 
when the best-fit line intensity was negative value, i.e. 3cr upper limit = max{best-fit line intensity, 0} 
-I- 3a statistical error range. This upper limit was shown in Eig. 3 (a) by the red line. We found line-like 
signatures at 0.600, 0.900, 1.275, 4.925 and 5.475 keV as shows in Eig. 3 (b), with 2 - 3a statistical 
significances. Since the energy of a DM emission line is not known a priori, it is necessary to take into 
account the number of independent line searches per one spectral fitting, i.e. the “trial factor”. Eor 
example, if we have a trial factor of 100 in a spectrum, the expectation to detect a line from random 
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Fig. 2. Exposure-time-weighted average of the 25 stacked XDB energy spectra from 2005 to 2013 and its best-fit model convolved with the corrected 
response. Sub-components of the model are represented by blue dashed (XDB in the Milky Way (1)+(2)+(4)), green dashed (CXB (3)) and orange dotted 
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Table 2. Spectral fitting resuits with the stacked energy spectra and the XDB model. 


Period 

Normi* * * * § 

kT2^ 

Norms* 

fcTst 

Ne* 

Mgt 


2005-2006 


u.z,z,_Q 03 

4 2~*"^'® 

^•^-0.4 

o.85i°:;« 

0.01J:§ 

o.olll 

2006-2007 

-7 + 2.5 
' -2.9 

^•^‘+-0.02 

9 -7+0.4 
' -0.4 

0.89t°:l° 

3.2lti 

6.5111 

2007-2008 



4 Q+4-2 
^•^-0.6 

0.79t°:l“ 

3.8l^-I 

0 1 

2008-2009 

20.0t?i®9 

fj '^-)+0.02 

U.Z,Z_o.o4 

4 <;+3.2 
^•^-0.5 

0.77t°;?? 

4.6l^;°s 

'^•^-0.0 

2009-2010 

14 q+11'3 

0 19+0.04 
U.l^_0.03 

6.8^®:? 

^•0/ -0.08 

2.911° 

0.5ll®5 

2010-2011 

42.2tti 

U.Z+_0.04 


0.69t°:“« 

1 .6lli 

o.ol°° 

2011-2012 

34.9t?4% 

rj '^'^+0.02 

U.Z,Z,_0_04 

c 0+3.6 

J.O-0.7 

0.641°;“ 

2.7111 

9 0+2.7 
■^‘°-2.0 

2012-2013 

36.8t«fo 

n 99+0.02 
U.z,z,_o 03 

Q q+3.8 
°-^-0.8 

o.55i°:;i 

3.3111 

4.0llj 


Period Norms* Pcxb^ •S'cxb'' OP x^^dof (dof) 


2005-2006 

0 

1 + 

0 0 

1 44+*^-*^^ 
^•^^-0.02 

7 

'•^-0.2 

0 . 411,1 

1.19 (590) 

2006-2007 

0 . 511,1 

1 49+0.02 

1 Q+0.2 
' '^-0.2 

0.611,1 

1.41 (434) 

2007-2008 

0.611,1 


9 O+0.2 

' •°-Q.2 

0.511,1 

1.22 (434) 

2008-2009 

0.611,1 

1 qq+0.02 

a.j:?_0.02 

-7 - 7 + 0.2 
'•'- 0.2 

0.411,1 

1.42 (434) 

2009-2010 

0.711° 

0 q 
d d 
+ 1 
00 

8 . 411,1 

0.211,1 

1.26 (434) 

2010-2011 

0.7111 

1 4 e+ 0-02 

a.h-o_q Q 2 

9 q+0.2 

'•^-0.2 

0.211,1 

1.44(434) 

2011-2012 

0.611,1 

i.JU_o.o2 

0 9+0.2 

O.Z-0.2 

4 9+0.3 
^•^- 0.3 

1.19 (434) 

2012-2013 

0.811,1 

1 or+0.02 
i.OJ_0_02 

9 O+0.2 
'•*^-0.2 

r O+ 0.4 
J-J- 0.4 

0.99 (434) 


Notes. 

* The emission measure of the optically-thin thermal CIE plasma integrated over the line of sight (the APEC model 
normalization): (l/27r) f Uenuds in unit of 10^“* cm“® sr“^, where rie and nn are the electron and the hydrogen density (cm“®). 

t The temperature of the optically-thin thermal CIE plasma in unit of ke V. 

J The abundances of Ne or Mg in unit of the Solar-neighbor values given in Anders & Grevesse (1989). 

§ The photon index of the power-law model for the CXB component. 

II The surface brightness of the CXB component (the power-law model normalization) in unit of photons cm“^ s“^ sr“^ keV“^ at 
1 keV (the photon index is fixed at 1.4). 

# The intensity of neutral oxygen (Oi, centroid: 0.525 keV) in unit of LU (photons cm“^ s“^ sr“^). 
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noise with an intensity exeeeding the 99 % signifieance is unity. This effeet is generally considered 
in the pulsation search from timing data (e.g. Woodet al. 1991), and is also sometimes called “Look- 
elsewhere Effect” (LEE; Gross & Vitells 2010). The trial factor is determined by the energy range and 
the energy resolution. This is not analytically known, since the energy spread function is relatively 
broad and since the energy resolution is energy-dependent. We therefore conducted Monte-Carlo 
simulations. We first generated a simulation spectrum assuming the best-fit XDB model without an 
additional emission line. We then determined the maximum significance for a dummy emission line 
in each simulated spectrum. This procedure was iterated 4000 times, and the cumulative occurrence 
distribution of the 4000 maximum significances was plotted in Eig. 4. Erom this figure, we found that 
significances > 3 .Oct appeared 635 times in the 4000 simulations (i.e. 16 %; the one-sided tail of p- 
value for Icr). This means that LEE-uncorrected 3. Oct corresponds to LEE-corrected la. Similarly, we 
found that LEE-uncorrected 4.la corresponded to LEE-corrected 3a. As a result, line-like signatures 
found in this search should be regarded as statistical fluctuations. 



Energy [keV] 



Energy [keV] 

Fig. 3. (a) Expected DM line intensity, its 1 <t statistical error range (blue crosses) and its So- upper limit (red line). The unit of line intensity is defined as “LU”, 
which is equal to photons cm”^ s“^ sr“^. (b) Five possible signatures and their statistical significances. Note that these significances are not and should 

be corrected by considering the LEE. 


We also evaluated the systematic uncertainty in this line search. Three main causes were 
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LEE-uncorrected Significance [a] 


Fig. 4. Probability distribution of mis-detection by statistical fluctuation with the 4000 simulations. The abscissa is the maximum LEE-uncorrected 
significance for a dummy emission line by random noise in a simulated spectrum, and the ordinate shows the number of occurrence in the 4000 simulated 
spectra, which is converted into the probability of occurrence (left scale) and the sigma value (right scale; LEE-corrected significance). The top 0.14, 2.3 and 
16 % of all (corresponding to the one-sided tail of p-value for 3, 2 and lo-) are distributed over 4.1,3.6 and S.Ocr, respectively. 


investigated: 

1. XDB model uneertainty. 

2. Instrumental response uneertainty. 

3. NXB eontribution uncertainty. 

The XDB model especially below 1 keV contains multiple emission lines from the hot plasma, i.e. the 
SWCX+LHB blend, the MWH and the UHTP, and their element abundances are difficult to determine. 
We varied their element abundances from 1/2 to 2 x the Solar-neighbor values and checked the effects 
on the line search results. The instrumental responses between 1.5 and 3.5 keV were corrected by the 
Crab calibration. The statistical uncertainty of this correction were also considered to give systematic 
errors. For the NXB reconstruction, the errors of the five instrumental lines were evaluated as the 
standard deviation of those obtained by the night-Earth data corresponds to the 187 blank sky data 
used in this analysis. 
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6 Discussion 


No significant signature of radiatively decaying DM was found in the 0.5 - 7 keV range in this search. 
This result is consistent with the stacked spectral analysis of the Perseus cluster by using Suzaku XIS 
observational data (Tamura et al. 2015). 

The DM line intensity / is expressed by the DM column density S'dm> the decay rate T and 
the mass moM: 

‘S'dmT 


1 = 


( 1 ) 


dvrrrtDM 

We are able to convert the DM line intensity normalized by its column density into its decay rate 
divided by its mass. We obtained the upper limit on the ratio of the DM decay rate and its mass 
as shown in Fig. 5 in comparison with previous works (Boyarsky et al. 2012; Bulbul et al. 2014; 
Boyarsky et al. 2014a). The exposure-time-weighted average of DM column density in this work was 
63 Mq pc“^ from the NFW model (Sofue 2012). The LEE-uncorrected upper limit was also plotted in 
addition to the EEE-corrected one because the previous works have not adapted the EEE-correction, 
and for the 3.5 keV upper limit, it does not need to be considered. Possible systematic uncertainties 
were also included. 

Assuming sterile neutrinos to make up DM, we converted this observational result into con¬ 
straints on their mass and mixing angle with the active neutrinos as the same manner in Abazajian 
et al. (2001a). The decay rate of sterile neutrinos is written as 


r = 


10247r2 


m®sin2 20 = 1.4x10-32 


ms \^/sin 2 26*' 


1 keV 




10 


-10 


( 2 ) 


where a is the fine-structure constant, Gp is the Eermi constant and 0 is a sterile neutrino mixing 
angle (Pal & Wolfenstein 1982). Prom Eq. 1 and 2, the line intensity is 

ms Y (sin^2e\ f / Sbm ^ 


1 = 1.3x10-3 LU 


(3) 


.1 keV; V 10”^° / V 1 / 

where /s is the fraction of Ug in DM. Prom Eq. 3, we obtained the constraint on their mass and 
mixing angle as shown in Pig. 6. The cyan and yellow shaded regions have been already rejected 
by several studies, such as the non-resonant (cyan shaded region; Boyarsky et al. 2009a), the reso¬ 
nant production with the maximal lepton asymmetry attainable in the z/MSM (yellow shaded region; 
Shaposhnikov 2008; Paine & Shaposhnikov 2008). The Tremaine-Gunn phase-space density consid¬ 
erations (Boyarsky et al. 2009a) and the Eyman-a analysis (Boyarsky et al. 2009b; Boyarsky et al. 
2009c) ruled out the region below 1 keV. We have added an exclude region for sin^ 26* in 3 < mg < 14 
keV. Horiuchi et al. (2014) also claimed a tight upper limit by using Chandra ACIS observational data 
of M31 with exposure time of ~ 400 ksec. We did not include it in our plot because their confidence 
level definition with fit statistics was not clear. 
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Energy [keV] 


- LEE-corrected 3o statistical + systematic upper limit 

- LEE-corrected 3o statistical upper limit 

- LEE-uncorrected 3 a statistical upper limit 


Fig. 5. Upper limit on the DM line intensity normalized by its column density corresponding to the ratio of DM decay rate and its mass. The LEE-uncorrected, 
-corrected 3a- and the LEE-corrected 3a statistical + systematic upper limits are indicated by the black dotted, dashed and solid lines, respectively. The 
typical 3a upper limit by previous works (Boyarsky et al. 2012) is expressed by the red dotted line. The parameters of the possible 3.5 keV line (Bulbul et al. 

2014; Boyarsky et al. 2014a) are plotted by the blue and magenta crosses. 


7 Conclusion and future prospects 

In this paper, we searched for the signature of an X-ray emission line from DM associated with the 
Milky Way halo by using a set of 187 Suzaku XIS archival data sets of the XDB from 2005 to 2013. 
After improving reproducibility of the instrumental responses and the instrumental line contributions, 
we searched for a non-baryonic emission line in the stacked XDB spectra by spectral fitting with [(the 
corrected response by using the stacked Crab spectra) x (the best-fit XDB model -i- five-instrumental 
lines -I- Gaussian emission line model for a DM signature)]. Consequently, we did not detect a possible 
DM signature including the line at 3.5 keV reported by previous studies (e.g. Bulbul et al. 2014), and 
determined the upper limit on the emission line intensity taking into account the LEE. We tightened 
the constraints on the ratio of DM decay rate and its mass in the 0.5 - 7 keV range and the parameters 
of sterile neutrinos. 

In the future, progressive instruments such as X-ray micro-calorimeters with energy resolution 
for diffuse X-ray emission of an order of eV and the large-EoV telescopes will be introduced to X- 
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fWs [keV] 

- LEE-corrected 3a statistical + systematic upper limit 

- LEE-corrected 3a statistical upper limit 

- LEE-uncorrected 3a statistical upper limit 

Fig. 6. Constraints on the sterile neutrino mass ms and mixing angle sin^ 29 by this and previous works. The definition of lines and marks are the same as 
in Fig. 5. The cyan and yellow shaded regions are excluded by the production theories of sterile neutrinos in the i^MSM (see text for details). 


ray observational satellites, and more sensitive DM searehes will be performed. Among the more 
hopeful instruments of near future missions, we are foeussed on the Soft X-ray Spectrometer (SXS) 
of the ASTRO-H satellite (Takahashi et al. 2010; Mitsuda et al. 2014) and the extended ROentgen 
Survey with an Imaging Telescope Array telescope with the PN-CCD camera module (eROSITA) of 
the Spektrum-Roentgen-Gamma (SRG) satellite (Predehl et al. 2014; Meidinger et al. 2014). The 
ASTRO-H SXS is an X-ray micro-calorimeter with doped semiconductor thermistors and will have 
the highest energy resolution, other than for grating instruments for point sources, though its grasp 
is lower than that of the existing X-ray observatories. On the other hand, the SRG eROSITA has the 
largest grasp and an all sky survey plan which is suitable for deeper analysis of the XDB, though its 
energy resolution is more modest. Especially in the ASTRO-H SXS with the high line identification 
ability by the high energy resolution, it is suitable for a weak line search with “dense” targets such as 
clusters of galaxies and nearby galaxies although their background plasma emission are strong. The 
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future X-ray observations will give a tighter eonstraint on DM eonditions. 


Acknowledgments 

This work is partly supported by JSPS KAKENHI Grant Numbers 12J10673, 25247028 and 26220703. We thank Dr. K. Sakai for valuable suggestions. 


References 

Abazajian, K., Fuller, G. M., & Tucker, W. H. 2001a, ApJ, 562, 593 
Abazajian, K., Fuller, G. M., & Patel, M. 2001b, Phys. Rev. D, 64, 023501 
Anders, E., & Grevesse, N. 1989, Geochim. Cosmochim. Acta, 53, 197 
Asaka, T., & Shaposhnikov, M. 2005, Physics Letters B, 620, 17 
Asaka, T., Blanchet, S., & Shaposhnikov, M. 2005, Physics Letters B, 631, 151 
Boyarsky, A., Nevalainen, J., & Ruchayskiy, O. 2007, A&A, 471, 51 

Boyarsky, A., Ruchayskiy, O., & lakubovskyi, D. 2009a, Journal of Cosmology and Astroparticle Physics, 3, 
5 

Boyarsky, A., Lesgourgues, J., Ruchayskiy, O., & Viel, M. 2009b, Journal of Cosmology and Astroparticle 
Physics, 5, 12 

Boyarsky, A., Lesgourgues, J., Ruchayskiy, O., & Viel, M. 2009c, Physical Review Letters, 102, 201304 
Boyarsky, A., Neronov, A., Ruchayskiy, O., & Tkachev, I. 2010, Physical Review Letters, 104, 191301 
Boyarsky, A., lakubovskyi, D., & Ruchayskiy, O. 2012, Physics of the Dark Universe, 1, 136 
Boyarsky, A., Ruchayskiy, O., lakubovskyi, D., & Lranse, J. 2014a, Physical Review Letters, 113, 251301 
Boyarsky, A., Lranse, J., lakubovskyi, D., & Ruchayskiy, O. 2014b, arxiv:1408.2503vl 
Bulbul, E., Markevitch, M., Loster, A., Smith, R. K., Loewenstein, M., & Randall, S. W. 2014, ApJ, 789, 13 
Carlson, E., Jeltema, T, & Profumo, S. 2015, Journal of Cosmology and Astroparticle Physics, 2, 9 
Cox, D. P, & Anderson, P R. 1982, ApJ, 253, 268 

Cox, D. P . 1998, in Lecture Notes in Physics, Berlin Springer Verlag, Vol. 506, lAU Colloq. 166: The Local 
Bubble and Beyond, ed. D. Breitschwerdt, M. J. Lreyberg, & J. Truemper, 121-131 
Cravens, T. E. 2000, ApJL, 532, L153 

Dodelson, S., & Widrow, L. M. 1994, Physical Review Letters, 72, 17 
Dolgov, A. D., & Hansen, S. H. 2002, Astroparticle Physics, 16, 339 
Poster, A. R., Ji, L., Smith, R. K., & Brickhouse, N. S. 2012, ApJ, 756, 128 
Pujimoto, R., et al. 2007, PASJ, 59, 133 
Pukuda, Y., et al. 1998, Physical Review Letters, 81, 1562 

18 


Gilmore, G., Wilkinson, M. I., Wyse, R. F. G., Kleyna, J. T., Koch, A., Evans, N. W., & Grebel, E. K. 2007, 
ApJ, 663, 948 

Goerdt, T., Moore, B., Read, J. I., Stadel, J., & Zemp, M. 2006, MNRAS, 368, 1073 
Gross, E., & Vitells, O. 2010, European Physical Journal C, 70, 525 

Hagihara, T., Yamasaki, N. Y, Mitsuda, K., Takei, Y, Sakai, K., Yao, Y, Wang, Q. D., & McCammon, D. 
2011, PASJ, 63,889 

Horiuchi, S., Humphrey, P J., Onorhe, J., Ahazajian, K. N., Kaplinghat, M., & Garrison-Kimmel, S. 2014, 
Phys. Rev. D, 89, 025017 
Ishisaki, Y, et al. 2007, PASJ, 59, 113 
Jeltema, T. and Profumo, S. 2015 MNRAS, 450, 2143 

Kalherla, P M. W., Burton, W. B., Hartmann, D., Arnal, E. M., Bajaja, E., Morras, R., & Pdppel, W. G. E. 

2005, A&A, 440, 775 

Koutroumpa, D., Eallement, R., Kharchenko, V., Dalgarno, A., Pepino, R., Izmodenov, V., & Quemerais, E. 

2006, A&A, 460, 289 
Koyama, K., et al. 2007, PASJ, 59, 23 

Kushino, A., Ishisaki, Y, Morita, U., Yamasaki, N. Y, Ishida, M., Ohashi, T., & Ueda, Y. 2002, PASJ, 54, 327 
Kusenko, A. . 2006, Physical Review Eetters, 97, 241301 

Eaine, M., & Shaposhnikov, M. 2008,Journal of Cosmology and Astroparticle Physics, 6, 31 
Eovell, M. R., Erenk, C. S., Eke, V. R., Jenkins, A., Gao, E., & Theuns, T. 2014, MNRAS, 439, 300 
Malyshev, D., Nerenov, A., & Eckert, D. 2014, Phys. Rev. D, 90, 103506 

Masui, K., Mitsuda, K., Yamasaki, N. Y, Takei, Y, Kimura, S., Yoshino, T., & McCammon, D. 2009, PASJ, 
61,115 

McCammon, D., & Sanders, W. T. 1990, ARA&A, 28, 657 

Meidinger, N., et al. 2014, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, 
Vol. 9144, Society of Photo-Optical Instrumentation Engineers (SPIE) 

Mitsuda, K., et al. 2007, PASJ, 59, 1 

Mitsuda, K., et al. 2014, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 

9144, Society of Photo-Optical Instrumentation Engineers (SPIE) 

Moretti, A., Campana, S., Eazzati, D., & Tagliaferri, G. 2003, ApJ, 588, 696 
Mushotzky, R. E, Cowie, E. L., Barger, A. J., & Arnaud, K. A. 2000, Nature, 404, 459 
Pal, P. B., & Wolfenstein, E. 1982, Phys. Rev. D, 25, 766 
Petraki, K., & Kusenko, A. 2008, Phys. Rev. D, 77, 065014 

Predehl, R, et al. 2014, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 
9144, Society of Photo-Optical Instrumentation Engineers (SPIE) 


19 


Riemer-S0rensen, S. 2014, arXiv: 1405.7943 

Sakai, K., Yao, Y, Mitsuda, K., Yamasaki, N. Y, Wang, Q. D., Takei, Y, & McCammon, D. 2014, PASJ, 66, 
83 

Sekiya, N., Yamasaki, N. Y, Mitsuda, K., & Takei, Y 2014a, in Suzaku-MAXl 2014: Expanding the Frontiers 
of the X-ray Universe, ed. M. Ishida, R. Petre, & K. Mitsuda, 66 
Sekiya, N., Yamasaki, N. Y, Mitsuda, K., & Takei, Y 2014b, PASJ, 66, L3 
Serlemitsos, P J., et al. 2007, PASJ, 59, 9 

Shaposhnikov, M., & Tkachev, I. 2006, Physics Letters B, 639, 414 

Shaposhnikov, M. . 2008, Journal of High Energy Physics, 8, 8 

Shi, X., & Fuller, G. M. 1999, Physical Review Letters, 82, 2832 

Smith, R. K., Brickhouse, N. S., Liedahl, D. A., & Raymond, J. C. 2001, ApJL, 556, L91 

Sofue, Y 2012, PASJ, 64, 75 

Takahashi, T, et al. 2010, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, 
Vol. 7732, Society of Photo-Optical Instrumentation Engineers (SPIE) 

Tamura, T., lizuka, R., Maeda, Y, Mitsuda, K., & Yamasaki, N. Y. 2015, PASJ, 67, 23 
Tawa, N., et al. 2008, PASJ, 60, 11 

Wood, K.S., Norris, J.P., Hertz, P, Vaughan, B.A., Michelson, P.F., Mitsuda, K., Lewin, W.H.G., van Paradjis, 
J., Pennix, W., &van der Klis, M. 1991, ApJ, 379, 295 

Wyse, R. F. G., & Gilmore, G. 2008, in lAU Symposium, Vol. 244, lAU Symposium, ed. J. I. Davies & M. J. 
Disney, 44-52 

Yao, Y, Wang, Q. D., Hagihara, T, Mitsuda, K., McCammon, D., & Yamasaki, N. Y. 2009, ApJ, 690, 143 
Yoshino, T., Mitsuda, K., Yamasaki, N. Y, Takei, Y, Hagihara, T, Masui, K., Bauer, M., McCammon, D., 
Fujimoto, R., Wang, Q.D., & Yao, Y. 2009, PASJ, 61, 805 
Yoshitake, H., Sakai, K., Mitsuda, K., Yamasaki, N. Y, Takei, Y, & Yamamoto, R. 2013, PASJ, 65, 32 


20 


